Flow cytometry profiling of cellular immune response in COVID-19 infected, recovered and vaccinated individuals.

INTRODUCTION: SARS-CoV-2 has infected over 753 million individuals and caused more than 6.8 million deaths globally to date. COVID-19 disease severity has been associated with SARS-CoV-2 induced hyper inflammation and the immune correlation with its pathogenesis remains unclear. Acute viral infection is characterised by vigorous coordinated innate and adaptive activation, including an early cellular response that correlates well with the amplitude of virus specific humoral response. OBJECTIVE: The present study covers a wide spectrum of cellular immune response against COVID-19, irrespective of infection and vaccination. METHODS: We analysed immune status of (a) COVID-19 hospitalised patients including deceased and recovered patients, and compared with home isolated and non-infected healthy individuals, and (b) infected home isolated individuals with vaccinated individuals, using flow cytometry. We performed flow cytometry analysis of PBMCs to determine non-specific cell-mediated immune response. RESULTS: The immune response revealed extensive induction and activation of multiple immune lineages, including T and B cells, Th17 regulatory subsets and M1, M2 macrophages in deceased and hospitalised recovered patients, vaccinated and healthy individuals. Compromised immune cell expression was observed in deceased patients even in later stages, while expression was restored in hospitalised recovered patients and home isolated individuals. CONCLUSION: The findings associated with recovery and convalescence define a new signature of cellular immune response that persists in individuals with SARS-CoV-2 infection and vaccination. The findings will help in providing a better understanding of COVID-19 disease and will aid in developing better therapeutic strategies for treatment.

Evaluation of altered miRNA expression pattern to predict COVID-19 severity.

Outbreak of COVID-19 pandemic in December 2019 affected millions of people globally. After substantial research, several biomarkers for COVID-19 have been validated however no specific and reliable biomarker for the prognosis of patients with COVID-19 infection exists. Present study was designed to identify specific biomarkers to predict COVID-19 severity and tool for formulating treatment. A small cohort of subjects (n = 43) were enrolled and categorized in four study groups; Dead (n = 16), Severe (n = 10) and Moderate (n = 7) patients and healthy controls (n = 10). Small RNA sequencing was done on Illumina platform after isolation of microRNA from peripheral blood. Differential expression (DE) of miRNA (patients groups compared to control) revealed 118 down-regulated and 103 up-regulated known miRNAs with fold change (FC) expression ≥2 folds and p ≤ 0.05. DE miRNAs were then subjected to functional enrichment and network analysis. Bioinformatic analysis resulted in 31 miRNAs (24 Down-regulated; 7 up-regulated) significantly associated with COVID-19 having AUC>0.8 obtained from ROC curve. Seventeen out of 31 DE miRNAs have been linked to COVID-19 in previous studies. Three miRNAs, hsa-miR-147b-5p and hsa-miR-107 (down-regulated) and hsa-miR-1299 (up-regulated) showed significant unique DE in Dead patients. Another set of 4 miRNAs, hsa-miR-224-5p (down-regulated) and hsa-miR-4659b-3p, hsa-miR-495-3p and hsa-miR-335-3p were differentially up-regulated uniquely in Severe patients. Members of three miRNA families, hsa-miR-20, hsa-miR-32 and hsa-miR-548 were significantly down-regulated in all patients group in comparison to healthy controls. Thus a distinct miRNA expression profile was observed in Dead, Severe and Moderate COVID-19 patients. Present study suggests a panel of miRNAs which identified in COVID-19 patients and could be utilized as potential diagnostic biomarkers for predicting COVID-19 severity.

"Mitochondrial pathogenic mutations and metabolic alterations associated with COVID-19 disease severity".

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) caused global pandemic and drastically affected the humankind. Mitochondrial mutations have been found to be associated with several respiratory diseases. Missense mutation and pathogenic mitochondrial variants might unveil the potential involvement of the mitochondrial genome in coronavirus disease 2019 (COVID-19) pathogenesis. The present study aims to elucidate the role of mitochondrial DNA (mtDNA) mutations, mitochondrial haplogroup, and energy metabolism in disease severity. The study was performed on 58 subjects comprising COVID-19-positive (n = 42) and negative (n = 16) individuals. COVID-19-positive subjects were further categorized into severe deceased (SD), severe recovered (SR), moderate (Mo), and mild (Mi) patients, while COVID-19-negative subjects were healthy control (HC) for the study. High throughput next-generation sequencing was done to investigate mtDNA mutations and haplogroups. The computational approach was applied to study the effect of mtDNA mutations on protein secondary structure. Real time polymerase chain reaction was used for mtDNA copy number determination and mitochondrial function parameters were also analyzed. We found 15 mtDNA mutations in MT-ND5, MT-ND4, MT-ND2, and MT-COI genes uniquely associated with COVID-19 severity affecting the secondary structure of proteins in COVID-19-positive subjects. Haplogroup analysis suggests that mtDNA haplogroups M3d1a and W3a1b might be potentially associated with COVID-19 pathophysiology. The mitochondrial function parameters were significantly altered in severe patients (SD and SR; p < 0.05). No significant relationship was found between mtDNA mutations and oxidative stress markers (p > 0.05). The study highlights the importance of mitochondrial reprogramming in COVID-19 patients and may provide a feasible approach toward finding a path for therapeutic interventions to COVID-19 disease.

Development of robust, indigenous ELISA for detection of IgG antibodies against CoV-2 N and S proteins: mass screening.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) pandemic has adversely affected humankind and caused millions of deaths globally since January 2020. Robust and quick serological tests such as antibody detection assays for SARS-CoV-2 provide relevant information and aid in the process of vaccine development and diagnostics, as well as in sero-epidemiological monitoring of antibody response to the virus. The receptor-binding domain (RBD) of spike and nucleocapsid protein are specific targets for detecting SARS-CoV-2 antibodies. Here, we present the development of a stable spike (S) and nucleocapsid (N) protein-based ELISA antibody detection test "CoroSuchak," with 99% sensitivity, 98% specificity, cost-effective, and detection in a minimum time for serodiagnosis and mass screening of the population for antibodies against SARS-CoV-2. Blood samples were analyzed from 374 SARS-CoV-2 reverse transcription-polymerase chain reaction (RT-PCR) positive, 772 negative and asymptomatic, and 874 random groups of subjects. We found that the antibody titer was significantly higher (p < 0.0001) in infected and vaccinated group compared to the only vaccinated and only infected group. Using enzyme-linked immunosorbent assay (ELISA), we detected SARS-CoV-2 immunoglobulin G (IgG) antibodies in 118/123 (96%) infected individuals, 570/653 (87%) non-infected but vaccinated individuals, 231/237 (97%) individuals who were both infected and vaccinated, and 499/874 (57%) from randomly selected individuals from the first and second waves of the pandemic. Similarly in the third wave, 14/14 (100%) infected and 16/20 (80%) RT-PCR-negative but symptomatic subjects were detected. Thus, the highly sensitive and specific in-house developed ELISA antibody detection kit "CoroSuchak" is extremely useful to determine the seroprevalence of SARS-CoV-2 antibodies in the coronavirus-exposed population. KEY POINTS: •Indigenous kit using a combination of spike and nucleocapsid proteins and peptide sequences. •High sensitivity and specificity to detect variants. •Highly sensitive for mass screening.

Mitochondrial DNA mutations contribute to high altitude pulmonary edema via increased oxidative stress and metabolic reprogramming during hypobaric hypoxia.

High altitude pulmonary edema (HAPE) is experienced by non-acclimatized sea level individuals on exposure to high altitude hypoxic conditions. Available evidence suggests that genetic factors and perturbed mitochondrial redox status may play an important role in HAPE pathophysiology. However, the precise mechanism has not been fully understood. In the present study, sequencing of mitochondrial DNA (mtDNA) from HAPE subjects and acclimatized controls was performed to identify pathogenic mutations and to determine their role in HAPE. Hypobaric hypoxia induced oxidative stress and metabolic alterations were also assessed in HAPE subjects. mtDNA copy number, mitochondrial oxidative phosphorylation (mtOXPHOS) activity, mitochondrial biogenesis were measured to determine mitochondrial functions. The data revealed that the mutations in Complex I genes affects the secondary structure of protein in HAPE subjects. Further, increased oxidative stress during hypobaric hypoxia, reduced mitochondrial biogenesis and mtOXPHOS activity induced metabolic reprogramming appears to contribute to mitochondrial dysfunctions in HAPE individuals. Haplogroup analysis suggests that mtDNA haplogroup H2a2a1 has potential contribution in the pathobiology of HAPE in lowlanders. This study also suggests contribution of altered mitochondrial functions in HAPE susceptibility.

Mitochondrial DNA sequencing reveals association of variants and haplogroup M33a2'3 with High altitude pulmonary edema susceptibility in Indian male lowlanders.

High Altitude Pulmonary Edema (HAPE) is a threatening disorder caused due to acute exposure to high altitude above 3000 m. Apart from multiple factors involved, the genetic factors also play an important function in the pathogenesis of HAPE. This study aims to evaluate the role of mtDNA polymorphism and their association with haplogroup in understanding the etiology of HAPE. In this study, all the HAPE susceptible and acclimatized control subjects could be classified into nine haplogroups pertaining mostly to Macrohaplogroup M and U. The frequency of haplogroup M was significantly higher in HAPE susceptibles whereas the haplogroup M33a2'3 was found only in HAPE susceptibles. The variant G4491A and A4944G of MT-ND2, A14002G of MT-ND5, and C8562T of MT-ATP8, were definition site of haplogroup M33a2'3. The frequency of A10398G of MT-ND3, A8701G of MT-ATP6 and C14766T of MT-CYB genes were significantly higher in HAPE susceptibles. mtDNA copy number also plays a significant synergistic role in HAPE susceptibility. Our findings suggests that variants in MT-ND2 and MT-ND5 were predicted to confer decreased protein stability in HAPE susceptibles and in particular, highly conserved variants G4491A, A4944G and A14002G associated with haplogroup M33a2'3 may be the primary cause of susceptibility to HAPE in Indian male lowlanders.

Whole genome variant analysis in three ethnically diverse Indians.

India represents an amazing confluence of geographically, linguistically and socially disparate ethnic populations (Indian Genome Variation Consortium, J Genet 87:3-20, 2008). Understanding the genetic diversity of Indian population remains a daunting task. In this paper we present detailed analysis of genomic variations (high-depth coverage (~ 30×) using Illumina Hiseq 2000 platform) from three healthy Indian male individuals each belonging to three geographically delineated regions and linguistic phylum viz. high altitude region of Ladakh (Tibeto-Burman linguistic phylum), sub mountainous region of Kumaun (Indo-European linguistic phylum) and sea level region of Telangana (Dravidian linguistic phylum) for probing the extent of genetic diversity in our population. The sequencing analysis provided high quality data (~ 95% of the total reads aligned to the human reference genome for each sample) and very good alignment quality (> 80% of the filtered mapped reads had a quality score of 60). A total of 4.3, 3.7 and 4.3 million single nucleotide variations were identified in the genome of high altitude, sub mountainous and sea level respectively by comparing with human reference genome. Approximately 17.3, 18.2, 17.4% of the variants were unique in the three genomes. The study identified many novel variations in the three diverse genomes (132,970 in Ladakh, 112,317 in Kumaun and 128,881 in Telangana individual) and is an important resource for creating a baseline and a comprehensive catalogue of human genomic variation across the Indian as well as the Asian continent.

Polygenic study of endurance-associated genetic markers ACE I/D, ACTN3 Arg(R)577Ter(X), CKMM A/G NcoI and eNOS Glu(G)298Asp(T) in male Gorkha soldiers.

BACKGROUND: Gorkhas, a sub-mountainous population of the Himalayan region, are known for strength and bravery. In the present study when "Gorkha" is used without brackets, we are mentioning Gorkhas of Tibeto-Burman origin. Physical capability, strength and endurance are important components of fitness associated with genetic traits. The aim of this study was to examine the endurance potential of male Gorkha soldiers, based on endurance-related genetic markers ACE I/D, ACTN3 Arg (R)577Ter(X), CKMM A/G NcoI and eNOS Glu(G)298Asp(T). METHODS: Genotypic and allelic frequencies were determined in 374 male Gorkha soldiers (Tibeto-Burman). These frequencies were compared with frequencies obtained from Gorkha (Indo-Aryan), high-altitude natives (Tibeto-Burman) and Indian lowlanders (Indo-Aryan). "Total genotype score" (TGS) was calculated from accumulated combination of polymorphisms with maximum value "100" for theoretically "optimal" polygenic score. Probability of occurrence of "optimal" endurance profile was also determined. RESULTS: ACE II genotypic frequency was highest in Tamangs followed by Gurungs, Rais, Limbus and Magars. No statistical difference in genotypic and allelic frequency of ACTN3 Arg(R)577Ter(X) was noted within the groups. Rais showed the highest CKMM A allele frequency (0.908) compared to other Gorkha (Tibeto-Burman) groups. Limbus and Tamangs showed the highest eNOS G allele frequency (0.938 and 0.915, respectively) compared to that of other groups. Probability of male Gorkha soldiers possessing a theoretically optimal polygenic endurance profile for four candidate polymorphisms was ~3.35% (1 in 30). Four percent of the population of male Gorkha soldiers (15 in 374) exhibited an optimal TGS 100, and 16% exhibited TGS 87 for endurance compared to male Indian soldiers belonging to the lowland (Indo-Aryan) and Gorkha (Indo-Aryan) populations suggesting an overall more "favourable" polygenic profile in the male Gorkha soldier (Tibeto-Burman) population. CONCLUSIONS: This study presents evidence of higher frequency of endurance-associated genes in the Gorkhas implying thereby that such genetically endowed individuals from the population may be selected and trained for achieving excellence in endurance-related elite sports activities.

Genetic Diversity in Gorkhas: an Autosomal STR Study.

Genotyping of highly polymorphic autosomal short tandem repeat (STR) markers is a potent tool for elucidating genetic diversity. In the present study, fifteen autosomal STR markers were analyzed in unrelated healthy male Gorkha individuals (n = 98) serving in the Indian Army by using AmpFlSTR Identifiler Plus PCR Amplification Kit. In total, 138 alleles were observed with corresponding allele frequencies ranging from 0.005 to 0.469. The studied loci were in Hardy-Weinberg Equilibrium (HWE). Heterozygosity ranged from 0.602 to 0.867. The most polymorphic locus was Fibrinogen Alpha (FGA) chain which was also the most discriminating locus as expected. Neighbor Joining (NJ) tree and principal component analysis (PCA) plot clustered the Gorkhas with those of Nepal and other Tibeto-Burman population while lowlander Indian population formed separate cluster substantiating the closeness of the Gorkhas with the Tibeto-Burman linguistic phyla. Furthermore, the dataset of STR markers obtained in the study presents a valuable information source of STR DNA profiles from personnel for usage in disaster victim identification in military exigencies and adds to the Indian database of military soldiers and military hospital repository.